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Abstract 

We assessed the vulnerability to climate change of 156 rare plant species. The species were selected 
from the 1625 rare species in California to comprise eight rarity types, classified according to range 
size, population size, and habitat specificity. For each of the 156 species, we first assigned a climate 
change vulnerability score using life history attributes and species distribution models, as specified by 
the Climate Change Vulnerability Index (CCVI) of NatureServe. The resulting CCVI scores were 
extremely vulnerable (n = 2), highly vulnerable (n = 40), moderately vulnerable (n = 57), presumed 
stable (n = 32), increase likely (n = 16), and insufficient evidence (n = 9). Piperia yadonii Rand, 
Morgan & Ackerman and Mimulus purpureus A. L. Grant were the species scored as extremely 
vulnerable. There was no correlation of the CCVI scores with rarity type, suggesting that climate 
change vulnerability cannot be inferred by simple categorizations based on geographic range and 
habitat preference. Second, we conducted a follow-up species distribution model sensitivity analysis 
that showed that the modeling results were highly dependent upon both model algorithm and choice 
of predictor variables. However, 60 of the 156 species were predicted to have declines in climatic 
suitability, regardless of modeling technique. Third, as an independent assessment of vulnerability, we 
calculated the topographic complexity around known occurrences of each species. Species in 
topographically dissected landscapes may be less vulnerable to climate change because they can find 
suitable climates locally as climate changes. We found that topographic complexity varied 
substantially, even within a single CCVI score level, and therefore provides unique information on 
vulnerability. Our results can be used to guide monitoring, management, and conservation plans for 
rare plant species. 
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Climate change may negatively impact the flora 
of California, a biodiversity hotspot with over 2000 
endemic plant species (Myers et al. 2000). A 
changing climate may reduce and extirpate popu¬ 
lations (Pounds et al. 2006), cause species to 
migrate north and upslope (Parmesan 1996; Kelly 
and Goulden 2008; Loarie et al. 2009), advance 
flowering times, promote species invasion, increase 
disturbance (e.g., fire), and cause community 
reorganization (Walther et al. 2002; Burkett et al. 
2005). Several tools have been developed to identify 
which species and habitats are most imperiled by 
the negative impacts of climate change (Schnieder 
et al. 2007; Williams et al. 2008; Heller and 
Zavaleta 2009; EPA 2009; Byers and Norris 2011; 
Glick et al. 2011; Schlesinger et al. 2011), under the 
assumption that the world will continue to warm in 
the near term, even if emissions are immediately 
reduced (IPCC 2007). By identifying species or 
habitats most at risk from effects of climate change, 
conservation and management efforts can be 
targeted to reduce these impacts, such as by 
protecting existing habitat or through assisted 
migration (Hunter 2007; McLachlan et al. 2007). 

Vulnerability assessments provide a standard¬ 
ized method to assess sensitivity to climate 


change that is time-efficient, repeatable, and is 
directly comparable across species. Investigators 
can use vulnerability assessments to rank a list of 
species with regard to their relative expected 
sensitivity to shifts in climate (Gardali et al. 
2012). Most studies are focused upon individual 
species; however, recent work has also considered 
the vulnerability of landscapes (Klausmeyer et al. 
2011). For example, topographically complex 
landscapes may provide refugia or paths for 
movement to more suitable areas under changing 
temperature and moisture regimes (Hunter et al. 
1988; Beier and Brost 2010). Moreover, under¬ 
standing how regional and local processes inter¬ 
act to create spatial heterogeneity in climate may 
help predict the direction and rate of climate 
change (Ackerly et al. 2010). Further, assessments 
of geophysical diversity (i.e., the number of 
bedrock types) may be a useful alternative to 
species-level assessments, since high levels of 
geological diversity are often associated with 
habitat heterogeneity and species diversity (An¬ 
derson and Ferree 2010). 

Species-level vulnerability assessments are typ¬ 
ically based on intrinsic life history traits, species 
distribution models (SDMs), or both. The trait- 
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based approach identifies and scores species 
attributes relevant to avoiding or tolerating 
climate change, such as movement ability (i.e., 
dispersal rates) or sensitivity to changes in 
temperature or moisture. The sum of these scores 
represents the species’ overall vulnerability to 
climate change. Trait-based indices were largely 
developed with animals as a primary focus. For 
example, the International Union for Conserva¬ 
tion of Nature (IUCN) analyzed the life history, 
ecology, behavior, physiology, and genetic make¬ 
up of “red list” animal species to assess a species 
vulnerability to climate change (Foden et al. 2009). 

A second set of tools used to assess vulnerabil¬ 
ity to climate change is SDMs (Pearson and 
Dawson 2003; Loarie et al. 2008; Stralberg et al. 
2009). Typically, point occurrence data for a 
species are used to create a statistical model of 
climatic suitability using historical (often 30-year- 
mean) climate data. This model is then used to 
predict the species’ contemporary range, based on 
a grid of historical climate, and the species’ future 
range, based on a grid of predicted future climate. 
Finally, the change in predicted range size and the 
amount of range overlap is calculated. Species 
with large range reductions and/or low range 
overlap are considered to be more vulnerable than 
species with small range reductions and/or high 
range overlap. However, there are a large number 
of modeling techniques used to describe climatic 
suitability, and differences in model algorithms 
and assumptions can greatly influence the quality 
of model predictions (Araujo and New 2007). 

A leading example of combining species traits 
and SDMs into a single vulnerability analysis 
comes from NatureServe (Arlington, VA), a non¬ 
profit organization whose mission is to provide 
the scientific basis for effective conservation 
action. NatureServe developed the Climate 
Change Vulnerability Index (CCVI) to serve as 
a standardized methodology for assessing vulner¬ 
ability to climate change at the species level 
(Young et al. 2012). The CCVI consists of a 
Microsoft Excel document with four main 
sections: Section A—direct exposure to changing 
temperature and precipitation; Section B—indi¬ 
rect exposure to climate change, including sea 
level rise, natural and human barriers, and land 
impacts from climate mitigation; Section C— 
sensitivity factors (hereafter referred to as “life 
history traits”); and Section D—modeled re¬ 
sponse to climate change. 

Our goal was to assess the vulnerability of 
California rare plant species to climate change 
and to evaluate the application of the Nature¬ 
Serve CCVI method to rare plants. Rarity is a 
major feature of California’s botanical heritage. 
The California Native Plant Society (CNPS) 
Rare Plant Program, which works in coordina¬ 
tion with the California Department of Fish and 
Game’s (CDFG) Natural Diversity Database, 


recognizes 1625 plant taxa as rare or endangered, 
as of March 1, 2011. While many of the 1625 taxa 
are subspecies and varieties, and thus the words 
“taxa” and “taxon” are appropriate than “spe¬ 
cies,” we use “species” throughout the text for 
simplicity. There are 26 California Rare Plant 
Rank 1A—presumed extinct in California; 1132 
Rank IB—rare or endangered in California and 
elsewhere; and 492 Rank 2—rare or endangered 
in California, but more common elsewhere 
(CNPS 2001). These rare species may have 
narrow ranges, small population sizes, or narrow 
habitat preferences (or all of the above) for 
natural or anthropogenic reasons (CNPS 2001). 
The potential impacts of climate change were not 
a factor considered by CNPS when assigning rare 
plant ranks; thus, which of the 1650 species will 
be most vulnerable to climate change has been 
largely uninvestigated. 

Climate-only SDMs suggest California plants 
may be in trouble: 66% will experience 80% 
reductions in range size within a century (Loarie 
et al. 2008). In addition, plants may be unable to 
adjust their ranges fast enough to spatially track 
shifting climates (Loarie et al. 2009). It is possible 
that rare species may be even more sensitive to 
climate change, given their limited geographic 
ranges and small population sizes. Or, perhaps 
they will be less sensitive to climate change, given 
their specialized ecologies. These attributes also 
make it more difficult to accurately model rare 
species than common species. Given our conser¬ 
vation concerns, it is our hope to create a 
meaningful vulnerability ranking for rare species 
and to identify which spatial and life history 
factors contribute most to that vulnerability. 

Due to the large number of rare plants in 
California, we sought to determine whether the 
level of climate change vulnerability could be 
inferred for certain groups of rare plants based on 
rarity type, life history traits, or biogeographic 
affinity. Our work can be divided into three 
complementary parts. First, for a subset of the 
1650 rare plants (n = 156), we compiled a set of life 
history attributes and created distribution models 
to rank vulnerability as specified by the CCVI of 
NatureServe. Second, we conducted a SDM 
sensitivity analysis to determine how choice of 
model algorithm and predictor variables influ¬ 
enced distribution model predictions of habitat 
suitability in future climates. Third, we calculated 
an index of vulnerability based on topographic 
complexity around known occurrences. 

Methods 
Species Selection 

To create a list of focal species that was 
representative of California rare plant species as a 
whole, we first classified each of the 1625 species 
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Table 1. Types of Rarity. The types of rarity, modified from Rabinowitz (1981), come from intersecting range 
size (small or large), population size (small or large), and habitat specificity (habitat specialist or generalist). For 
each category, the number of species in our sample and an example species is provided. While only seven of the 
eight groups are “rare,” species with relatively large ranges, large populations, and generalist habitat preferences 
among our sample are still relatively rare with respect to the average species in the flora. Thus, we used all eight 
groups for our selection. 



Large range 

Small 

range 

Population size 

Habitat generalist 

Habitat specialist 

Habitat generalist 

Habitat specialist 

Large populations 

Small populations 

n = 28 

California macrophylla 
n = 24 

Lilium parryi 

n = 25 

Streptanthus morrisonii 
n = 15 

Calochortus plummerae 

n = 24 

Mimulus purpureus 
n = 17 
Taraxacum 
calif or nicum 

n = 17 

Allium tuolumnense 
n = 6 
Monardella 
stebbinsii 


into one of the eight types or forms of rarity 
(Rabinowitz 1981) (Table 1). Following Rabino¬ 
witz’s (1981) definitions, only seven of the eight 
groups are “rare,” because the combination 
“large range, large population, and habitat 
generalist” is considered common. However, 
since our sample pool was made up of only rare 
species, even the species within this pool with 
large ranges, large populations, and generalist 
habitat preferences were rare relative to the 
average species in the flora. Thus, we sampled 
across all eight groups to obtain our subset of 156 
species (Appendix 1). Nomenclature for these 
species follows the CNPS Inventory (2001). 

We used information from the California 
Natural Diversity Data Base (CNDDB) to 
attribute each species with the three variables 
required for rarity type classification, as follows: 
(1) range size—total area of species range based 
on a minimum convex polygon encompassing 
mapped occurrences from the CNDDB; (2) 
population size—the median population number 
of individuals, extracted from the comment field 
of CNDDB; and (3) habitat specificity—sub¬ 
strate affinity, extracted from the habitat field of 
CNDDB. For range size and population size, 
species were designated as large or small based on 
their value relative to the median of the 
distribution of values. We then randomly selected 
species from each of the eight rarity types 
(Table 1). Our list of 156 species includes 139 
California Rare Plant Rank IBs, 13 Rank 2 s, 
and three Rank 3 s. While rank 3 s are not 
nominally rare and thus do not contribute to the 
1625 described above, we included three of these 
species to see if they had remarkably different 
vulnerability scores. 

Climate Change Vulnerability Index (CCVI) 

CCVI overview. The NatureServe CCVI (re¬ 
lease 2.01) assesses 24 climate change vulnerabil¬ 
ity risk factors, placed in four categories: direct 
exposure, indirect exposure, life history traits, 
and modeled response. To complete the CCVI, 
we collected information on the distribution. 


natural history, and conservation status of rare 
species from CDFG, CNPS, and NatureServe. 
We then conducted a literature review, mapped 
species distributions, and modeled responses to 
climate change. Biologists and botanists were 
consulted to fill data gaps as needed for 
particular species. Collected data and sources 
for each species are available online at www.dfg. 
ca.gov/biogeodata/. 

CCVI Section A: Direct exposure (two factors). 
Direct exposure was scored based on the 
percentage of the species’ range that falls into 
NatureServe’s recommended categories of pro¬ 
jected changes of temperature or moisture. The 
temperature change categories (decreasing in 
severity) were >3.1°C, 2.8-3.1°C, 2.5-2.7°C, 2.2- 
2.4°C, and <2.2°C; the moisture change categories 
(decreasing in severity) were <—0.119, —0.097- 
-0.119, -0.074— 0.096, -0.051—0.073, -0.028- 
—0.050, and >—0.028. Climate data (Fig. la-d) 
and projections for the year 2080 were derived by 
The Nature Conservancy and downloaded from 
their Climate Wizard (Model: Ensemble Average, 
emission scenario [ES] A2; www.climatewizard. 
org). Temperature change was the predicted 
change in annual temperature by 2080, calculated 
over the range of the species in California. 
Moisture change was the predicted net change in 
moisture based on the Hamon AET:PET Mois¬ 
ture Metric, calculated over the range of the 
species in California. Additional climate data was 
acquired for the modeled response (Section D) 
and SDM sensitivity analysis (II) from World- 
Clim, as described below. 

CCVI Section B: Indirect exposure (four 
factors). Indirect exposure evaluated landscape 
configuration factors that may affect the vulner¬ 
ability of a species to climate change: Exposure to 
sea level rise, distribution relative to natural 
barriers, distribution relative to anthropogenic 
barriers, and predicted impact of land use 
changes resulting specifically from human re¬ 
sponses to climate change. To evaluate these 
factors, we compared the distribution of the 
known occurrences of each species with a map of 
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Fig. 1. Map of future temperature change (a), future drying (b), historical temperature variation (c), and 
historical annual precipitation (d). 


predicted sea level rise (Strauss et al. 2012), 
topographic maps that depicted landscape fea¬ 
tures, and maps of proposed renewable energy 
(e.g., solar power stations, wind farms, geother¬ 
mal wells; California Department of Fish and 
Game 2011). For sea level rise, the percent of the 


species’ range occurring in areas subject to sea 
level rise were placed into five broad categories 
(decreasing in vulnerability): (1) >90%, (2) 50- 
90%, (3) 10—49, (4) <10, and (5) predicted 
increase in extent (intertidal species whose habitat 
may increase with sea level rise). For barriers. 
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Table 2. Life History Traits Used in Determining Climate Change Vulnerability Index (CCVI) in 
this Study. Assumptions and data sources are listed for each determining factor. 


Factor 


Assumption 


Data source 


Dispersal and movements 

Predicted sensitivity to 

temperature and moisture 
changes 


Dependence on a specific 
disturbance regime 

Restriction to uncommon 
geological features or 
derivatives 

Reliance on interspecific 
interactions 


Genetics 

Phenological response 


Lower dispersal ability leads to 
high vulnerability. 

Narrow historical climate exposure 
and special microclimatic 
preferences leads to high 
vulnerability. 


Dependence on a particular 
disturbance regime leads to 
high vulnerability. 

Habitat specialization leads to 
high vulnerability. 

Dependence on other species 
(facilitation, pollinators, 
and seed dispersers) leads to 
high vulnerability. 

Low genetic diversity leads to 
high vulnerability. 

Shorter bloom period leads to 
high vulnerability. 


Dispersal mechanisms from scientific 
literature and expert opinion. 

Historical temperature and 

precipitation variation from Climate 
Wizard. Physiological hydrological 
niche and physiological thermal 
niche from scientific literature and 
expert opinion. 

Adaptation and affinity for fire and 
flood-prone habitats from scientific 
literature and expert opinion. 

Substrate affinity from CNDDB. 
SSURGO soil data from the NRCS. 

Literature and expert opinion. 


Not scored for any species due to 
insufficient information. 

CNPS’s Rare Plant Inventory bloom- 
period database. 


each species was placed into one of four 
categories (decreasing in vulnerability): (1) Barri¬ 
ers completely surround the current distribution, 
(2) barriers border the current distribution 
incompletely but will likely impair distributional 
shits, (3) barriers border the current distribution 
incompletely but will be unlikely to impair 
distribution shifts, and (4) significant barriers 
do not exist. For renewable energy, each species 
was placed into one of four categories (decreasing 
in vulnerability): The likelihood that a species’ 
natural history or range may be (1) very likely to 
conflict with mitigation-related land-use changes, 
(2) likely to conflict with land-use changes, (3) 
likely to benefit from land-use changes, or (4) 
very likely to benefit from land-use changes. 

CCVI section C: life history traits (16 factors). 
The life history traits were grouped into the 
following categories: dispersal/movement, sensitiv¬ 
ity to temperature or moisture, disturbance-depen¬ 
dence, geologic restriction, interspecific interac¬ 
tions, genetic diversity, and phenology (Table 2). 
Each species was scored by assessing whether its life 
history traits would be expected to decrease, 
somewhat decrease, neutral, somewhat increase, 
increase, or greatly increase its vulnerability to 
climate change. If information was not available for 
a particular factor, it was scored as unknown. 
Guidance on how to rank each factor was provided 
by NatureServe. Ranks for ten factors in this group 
were required or the vulnerability index returned a 
score of “Insufficient Evidence.” 

CCVI Section D: Modeled response (two 
factors). We modeled the change in range size 


and range overlap of predicted future range with 
predicted current range using the Maxent algo¬ 
rithm, a statistical model that uses machine 
learning (Phillips et al. 2006; Elith and Leathwick 
2009). 

The spatial centroid of each CNDDB mapped 
occurrence record was used. The records of the 
CNDDB have been carefully curated by the staff 
of CDFG, yet uncertainty in the exact locations 
likely remains, contributing an unknown, but 
likely small, amount of variation to the model 
results presented here. Another limitation is that 
occurrences outside of CA are omitted; however, 
visual examination of the distributions of our 
focal species suggests that only —10% have 
distributions that may cross into Oregon or Baja 
California, Mexico. We used a dataset acquired 
from WorldClim (Hijmans et al. 2005) compris¬ 
ing four climate variables (annual temperature, 
annual precipitation, seasonality of temperature, 
and seasonality of precipitation) for current 
conditions (mean 1950-mean 2000) and for 
future conditions (called “2080,” but data are 
means for 2070-2100; Global Circulation Model 
[GCM] CGCM3.1, ES A IB) at 1 km 2 resolution. 
We fit a Maxent model for current conditions 
and used the resulting model to predict climatic 
suitability, ranging from 0-1, for both current 
and future conditions; background points were 
selected randomly. Maxent was run using the 
default “auto features” mode, allowing the use of 
linear, quadratic, product, threshold, and hinge 
features. These continuous surfaces were then 
converted to binary (suitable/unsuitable) using a 
threshold determined as the value that maximizes 
the kappa, a statistical measure of the agreement 





198 


MADRONO 


[Vol. 60 


Table 3. Climate Change Vulnerability Index Score Descriptions. 


Index scores 


Descriptions 


Extremely vulnerable (EV) 
Highly vulnerable (HV) 
Moderately vulnerable (MV) 
Presumed stable (PS) 

Increase likely (IL) 
Insufficient evidence (IE) 


Abundance and/or range extent within geographical area assessed 
extremely likely to substantially decrease or disappear by 2050. 

Abundance and/or range extent within geographical area assessed likely to 
decrease significantly by 2050. 

Abundance and/or range extent within geographical area assessed likely to 
decrease by 2050. 

Available evidence does not suggest that abundance and/or range extent 
within the geographical area assessed will change (increase/decrease) 
substantially by 2050. Actual range boundaries may change. 

Available evidence suggests that abundance and/or range extent within 
geographical area assessed is likely to increase by 2050. 

Available information about a species’ vulnerability is inadequate to 
calculate an index score. 


between predictions and observations, and the 
AUC values were calculated (Cohen 1960; 
Jimenez-Valverde 2011). Change in range size 
was calculated as the total area predicted as 
suitable in the future (t 2 ) minus the total area 
predicted as suitable in the present (ti), divided by 
total area predicted as suitable at C. Range 
overlap was calculated as the total area predicted 
to be suitable at ti and t 2 , divided by the total 
area predicted as suitable at tp The geographic 
extent of all models was California. This may 
overestimate range loss when a species’ new range 
is predicted to be outside of California. 

CCVI risk factor score. The natural history and 
distributional information for each species was 
entered into the CCVI Excel calculator to obtain 
scores for each species. The output was one of six 
vulnerability scores: extremely vulnerable (EV), 
highly vulnerably (HV), moderately vulnerable 
(MV), presumed stable (PS), increase likely (IL), 
and insufficient evidence (IE; Table 3). All 
vulnerability index scores were calculated with 
and without modeled response to climate change 
(Section D). 

Statistical analysis of CCVI predictors. To 
identify which of the factors were most strongly 
associated with the resulting vulnerability scores, 
we evaluated the distribution of risk factor scores 
against the distribution of CCVI scores across all 
of our species. Factors that were frequently scored 
as increasing or decreasing vulnerability should 
show up as significant predictors of the distribu¬ 
tion of CCVI scores, while factors that were 
infrequently scored as increasing or decreasing 
vulnerability should not. We converted the CCVI 
scores and factor ranks to their numeric equiva¬ 
lents, then regressed the CCVI score vs. each of 
the factors from these sections. We did not test for 
a relationship of CCVI with dietary versatility, 
genetic variation, or genetic bottlenecks, because 
no species were scored for those factors. Tests 
varied in the number of species included because 
we did not have information on all life history 


traits for all species (i.e., when the factors were 
marked as unknown). We also compared the 
CCVI score with California Rare Plant Rank and 
rarity type using two one-way ANOVAs. Finally, 
we tested if range size change or range overlap 
was related to California Rare Plant Rank and 
rarity type using four one-way ANOVAs. 

SDM Sensitivity Analysis 

For each of the 156 species, we ran 22 
additional SDMs to estimate the sensitivity of 
range predictions to modeling algorithms and 
choice of predictor variables. Our 23 models (all 
run in Maxent unless stated otherwise) were: 1) 
19 climate variables (bioclim); 2) four climate 
variables (described above); 3-14) four climate 
variables, with different GCM*ES combinations 
(GCMs included BCCR-BCM2.0, CSIRO- 
Mk3.0, INM-CM3.0, and MIROC3.2 [medres]; 
ESs included A2, AB, and Bl); 15) 19 climate 
variables with soil type; 16) 19 climate variables 
with soil properties (pH, organic matter, and 
clay); 17) four climate variables with soil type; 18) 
four climate variables with soil properties (pH, 
organic matter, and clay); 19) four climate 
variables, with a customized geographic extent 
for each species; 20) four climate variables with 
an equal number of presences and psuedo- 
absences; 21) four climate variables, with Ran¬ 
dom Forest; 22) four climate variables, with 
random forest, and with a customized geographic 
extent for each species; 23) four climate variables, 
with the boosted regression tree model. Soil type 
data came from the Geologic Map of California 
(Jennings et al. 2010), which we simplified into 
seven “soil types” (gabbro, granite, limestone, 
sandstone, serpentine, shale, and volcanic) and 
rasterized to 1 km 2 resolution. Soil property data 
were obtained from the Natural Resources 
Conservation Service (NRCS) soil survey geo¬ 
graphic database (SSURGO) and rasterized to 
1 km 2 . The customized geographic extents were 
determined by intersecting the point occurrence 
data with the Jepson Ecoregions (Hickman 1993); 
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A. Eriogonum twisselmannii 


B.Limosella australis 


A 



Fig. 2. Illustration of topographic complexity, (a) Eriogonum twisselmannii had an anomaly score of —0.61 and a 
topographic complexity (standard deviation of elevation) of 27.3. (b) Limosella australis had an anomaly score of 
—0.64 and a topographic complexity score of 0.51. 


for each species, the predictor raster surfaces were 
cropped to the shape of the ecoregions that 
contained the species. 

To compare the model predictions, we calculated 
an “anomaly score” for each species and each model. 
For a given species, at each known occurrence, an 
anomaly value was calculated as the predicted 
suitability in the present (ti) subtracted from the 
predicted suitability in the future (t 2 ). The mean of 
the anomaly values across all occurrences was the 
“anomaly score” for each species. A negative 
anomaly score meant that suitability in the future 
was predicted to be lower than current suitability. 

There are two main advantages of using an 
anomaly score to compare model predictions for 
our SDM sensitivity analysis, rather than using 
the conventional metrics of change in range size 
or range overlap. First, the anomaly score can be 
derived without converting continuous predicted 
surfaces to binary, the latter of which requires 
determining a threshold above or below which a 
particular location is considered suitable or 
unsuitable. For rare species, it is especially 
difficult to choose a meaningful threshold, given 
the limited number of point occurrences available 
for evaluation of the prediction errors made at 
various threshold values (Hijmans 2012). Second, 
anomaly scores are based on climate suitability 
change only at known occurrences, while range 
size and overlap consider the entirety of Califor¬ 
nia. Most rare plant species occupy only a small 
portion of their range where specific habitat 
requirements are met, and managers are most 
concerned with how climate suitability will 
change where the species actually occurs. Evalu¬ 
ating change in suitability over the entire range of 
the species requires assumptions about habitat 
occupancy and movement that may not be met 
by many rare plant species. Furthermore, using 
the anomaly score meant that predictions were 
only necessary for known occurrences, making 
our sensitivity analysis computationally feasible. 

We asked if the anomaly score was significant¬ 
ly related to the model type, within species, using 
a two-way ANOVA. Tukey’s HSD tests were 


used for post-hoc means separation by model 
type. We also examined the relationship of 
anomaly score and CCVI score using linear 
regression. Finally, we tested if the median 
anomaly score of each species was related to 
California Rare Plant Rank or rarity type using 
two one-way ANOVAs. 

Topographic Complexity Analysis 

As an index of “topographic complexity,” we 
calculated the standard deviation of elevation 
(resolution = 30 X 30 m) within 100 m of each 
occurrence, and took the mean of those values 
per species. Topographic complexity measured at 
1000 m was highly correlated with topographic 
complexity measured at 100 m (r = 0.94, P < 
0.001), so just the 100 m buffer was used. A 
species in a topographically complex landscape 
was considered less vulnerable than a species in a 
topographically homogeneous landscape (Fig. 2). 
We asked if the topographic complexity score 
was significantly related to the CCVI score using 
linear regression. We also tested if the topo¬ 
graphic complexity score was related to Califor¬ 
nia Rare Plant Rank or rarity type using two one¬ 
way ANOVAs. 

Results 

Climate Change Vulnerability Index 

Of the 156 species assessed, 99 were determined 
to be vulnerable (extremely vulnerable, highly 
vulnerable, or moderately vulnerable) to climate 
change and 48 were determined to be stable or 
increasing (presumed stable or increase likely). 
The distribution of final scores was: extremely 
vulnerable (n = 2), highly vulnerable (n = 40), 
moderately vulnerable (n = 57), presumed stable 
(n = 32), increase likely (n = 16), and insufficient 
evidence (n = 9). Piperia yadonii Rand, Morgan 
& Ackerman and Mimuluspurpureus A. L. Grant 
were the species scored as extremely vulnerable. 
All assessment scores and species attribute data 
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Table 4. Top Five Most Vulnerable Species as Determined by Climate Change Vulnerability Index 
Score. For CCVI, IL = increase likely, PS = presumed stable, MV = moderately vulnerable, HV = highly 
vulnerable, EV = extremely vulnerable, and IE = insufficient evidence. CCVI (without D) is the CCVI score 
recalculated after excluding the species distribution model results. The global (G) and state (S) rankings are from 
the California Natural Diversity Data Base; G-rank is the global rank, reflecting global rarity, and S-rank is the 
state rank, reflecting state rarity. Both indices range from 1 to 5, with 1 being critically imperiled and 5 being secure. 
California Rare Plant Rank IB is rare or endangered in California and elsewhere and Rank 2 is rare or endangered 
in California, but more common elsewhere. Federal listing is the status of the species according to the Endangered 
Species Act. Section D refers to the modeled response of the CCVI. CCVI index score definitions follow Table 3. 
*Limosella australis may not be native to California. 


Rank 

Species 

CCVI 

CCVI 

(without D) 

G-rank 

S-rank 

California 
Rare Plant 
Rank 

Federal 

listing 

Based on CCVI with Section D 

1 Piperia yadonii (Orchidaceae) 

EV 

HV 

5 

2 

2 

None 

2 

Mimulus purpureus (Phrymaceae) 

EV 

HV 

2 

2.2 

IB 

None 

3 

Calliandra eriophylla (Fabaceae) 

HV 

MV 

5 

2&3 

2 

None 

4 

Limosella australis* 

HV 

HV 

4&5 

2 

2 

None 

5 

(Scrophulariaceae) 

Taraxacum californicum 

HV 

MV 

2 

2 

IB 

Endangered 

Based 

1 

(Asteraceae) 

on CCVI without Section D 

Monolopia congdonii (Asteraceae) 

MV 

EV 

3 

3 

IB 

Endangered 

2 

Orcuttia viscida (Poaceae) 

HV 

EV 

1 

1 

IB 

Endangered 

3 

Pogogyne abramsii (Lamiaceae) 

MV 

EV 

1 

1 

IB 

Endangered 

4 

Symphyotrichum lentum 

HV 

EV 

2 

2 

IB 

None 

5 

(Asteraceae) 

Mimulus purpureus (Phrymaceae) 

EV 

HV 

2 

2.2 

IB 

None 


are reported in Appendix 1. The top five most 
vulnerable species, with and without Section D 
(“modeled response”), are listed in Table 4. 

For the 156 species assessed, the final CCVI 
score was significantly related to just one of the 
factors that were used to calculate it: anthropo¬ 
genic barriers (NatureServe factor B2b; r 2 = 0.09, 
P < 0.001, n = 147). Anthropogenic barriers 
were determined to limit the migration ability of 
99 of the 156 species, and the presence of 
anthropogenic barriers was significantly related 
to overall climate change vulnerability. Although 
not significantly related to overall climate change 
vulnerability, two factors were found to increase 
vulnerability for the majority of rare plants 
assessed: land use change from human response 
to climate change was found to increase vulner¬ 
ability for 80 species, and narrow temperature 
tolerance (“historical thermal niche”) was also 
found to increase vulnerability for 80 species. 

Overall climate change vulnerability was not 
significantly related to California Rare Plant 
Rank. For example, final CCVI scores for IB 
species were spread rather evenly as highly 
vulnerable (n = 35), moderately vulnerable (n 
= 49), presumed stable (n = 30), and increase 
likely (n = 14). In addition, overall climate 
change vulnerability was not related to rarity 
type. Similarly, neither change in range size nor 
range overlap were significantly related to Cali¬ 
fornia Rare Plant Rank or rarity type. 


SDM Sensitivity Analysis 

The anomaly scores statistically differed by 
species and by model type (species P < 0.001; 
model type P < 0.001), where the models with the 
most positive anomalies (greatest increase in 
suitability) were those made with Maxent that 
included soils information, and the models with 
the most negative anomalies (greatest reduction is 
suitability) were those made with Random 
Forest. 

The anomaly score and the CCVI score were 
significantly related (r 2 = 0.54, P < 0.001), where 
species scored as vulnerable also had low 
anomaly scores. This is not surprising, given the 
CCVI score included modeled response. In fact, 
the CCVI score calculated after excluding mod¬ 
eled response was not significantly related to the 
anomaly score. The anomaly score was not 
significantly related to California Rare Plant 
Rank or rarity type. 

Topographic Complexity Analysis 

Topographic complexity and the CCVI score 
were not significantly related. Topographic com¬ 
plexity was not significantly related to California 
Rare Plant Rank, but was significantly related to 
rarity type (P < 0.001), where habitat specialists 
occurred in locations with higher topographic 
complexity than habitat generalists. 
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Discussion 

We have three key results: (1) 2/3 of our focal 
species were scored as vulnerable to climate 
change, (2) SDM predictions were highly vari¬ 
able, and (3) topographic complexity may pro¬ 
vide complementary information on climate 
change vulnerability. 

Climate Change Vulnerability Index (CCVI) 

Ninety-nine of our 156 species (63%) were 
vulnerable to climate change (scored as moder¬ 
ately vulnerable or worse). We present a list of 
the top five most vulnerable species (Table 4), an 
annotated species list (Appendix 1), and all the 
information we used to make our determinations 
(CCVI Excel workbooks and species profiles 
available online; CDFG 2011). We were unable 
to elucidate strong relationships between species 
characters and vulnerability. No significant 
relationship was found between the CCVI 
vulnerability rank and California Rare Plant 
Rank, plant rarity type (Rabinowitz 1981), or 
any species life history trait considered. This 
suggests that direct exposure to climate change 
based on projected changes to future temperature 
and precipitation conditions within a species’ 
range was the strongest driver of vulnerability. 

One vulnerability factor from the Indirect 
Exposure Section (landscape configuration) was 
significantly (but weakly) related to the CCVI 
scores: anthropogenic barriers. For 99 of the 156 
species, the ability to migrate to track shifting 
climate will likely be impeded by man-made 
barriers. This may reflect the fact that many rare 
species are concentrated in coastal areas, where 
population density and associated fragmentation 
are extremely high (e.g., the San Francisco Bay 
Area and Los Angeles) (Stein et al. 2000). Coastal 
areas support “naturally rare” plant species that 
have evolved or survived over time in local 
refugia because of cool, aseasonal climates and 
a high level of soil heterogeneity. Other rare 
plants in coastal areas were once more common 
and have become “anthropogenically rare” due 
to high levels of development and habitat loss. 
Whether naturally or anthropogenically rare, if 
these species are unable to tolerate new climate 
conditions and cannot find refuge from novel, 
intolerable climates locally, the likelihood of 
dispersing to a more favorable, distant location 
is expected to be very low given man-made 
barriers. These species are prime candidates for 
assisted migration (McLachlan et al. 2007; 
Richardson et al. 2009; Vitt et al. 2010). 

The general inability to predict the CCVI 
vulnerability scores using the indirect exposure 
and life-history traits reflects the high importance 
of direct exposure (projected changes to temper¬ 
ature and moisture within a species’ range) in 


calculating the overall CCVI score. This is 
perhaps unsurprising, given that the indirect 
exposure and life-history traits are weighted by 
direct exposure in the calculation of the final 
CCVI score (Young et al. 2012). While this 
multiplicative approach makes it difficult to 
identify which risk factors are most important 
for a group of species, it is the appropriate 
approach for the CCVI, because a species with 
zero exposure to climate change should be 
considered invulnerable even if it has traits that 
make it sensitive to climate change, and vice 
versa. 

We found the CCVI scores were independent 
of both California Rare Plant Rank and rarity 
type. This suggests that the rarest species, such as 
the California Rare Plant Rank IB species or 
habitat specialists with small ranges and small 
population sizes, are not necessarily the most 
vulnerable to climate change. We also found no 
relationship of CCVI scores with other species 
attribute information, including plant life-cycle 
duration, plant growth form, and biogeographic 
affinity (results not shown). This is a desirable 
result, as it suggests that the CCVI scores contain 
novel information. This is perhaps the greatest 
strength of the CCVI: It represents a framework 
for thinking exclusively about climate change 
vulnerability. In developing the CCVI, a goal of 
NatureServe was to create an index that would be 
as independent as possible from existing rarity 
ranking indices. We suggest that the CCVI score 
can be viewed alongside of rarity scores (e.g., 
California Rare Plant Rank) to identify the most 
vulnerable and rarest species. 

SDM Sensitivity Analysis 

Our sensitivity analysis found that SDM 
predictions were highly variable. The range of 
anomaly scores per species was very large 
(Fig. 3), reflecting sensitivities to the choice of 
predictor variables and model algorithm (i.e., 
Maxent, Random Forest, and boosted regression 
trees). Further, for 90 of the 156 species, the 
direction of the anomaly (decreasing suitability or 
increasing suitability) was even variable (yellow 
bars of Fig. 3). One of the largest sources of 
variability in modeled response was the algo¬ 
rithm. In particular, Random Forest always 
produced the highest anomaly scores. Despite 
the high variability in anomaly scores, AUC 
values were uniformly high (mean = 0.988, 
median = 0.993, range = 0.921-1.00). These 
high AUC values suggest that all models are 
equally good, despite the fact that the models 
make radically different predictions about future 
climatic suitability. This is unsurprising, as AUC 
values tend to be inflated for species with narrow 
distributions, and therefore may not adequately 
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Fig. 3. Sensitivity analysis results, by species. The variation in anomaly score across 24 distribution models is 
represented with box plots, where open circles indicate statistical outliers. Dark grey boxes are those species whose 
distributions do not cross the horizontal dashed line. Light grey boxes are those species whose distributions do 
cross the horizontal dashed line. 


capture differences in model success for rare 
plants (Lobo et al. 2007; Gogol-Prokurat 2011). 

The variability in our modeled results is 
probably due in large part to working with rare 
species, which have few occurrences (the median 
number of occurrences for our species was 35) 
and many non-climatic factors that detennine 
their distributions. This reflects the “rare species 
modeling paradox,” which says that the species 
“most in need of predictive distribution modeling 
... are the most difficult to model” (Lomba et al. 
2010). In other words, while it is extremely 
important to identify the climatic tolerances of 
rare species, given our conservation concerns, it is 
extremely difficult to make meaningful models of 
the climatic ranges based on the limited available 
information. 

The variability of the modeled results is cause 
for concern because our CCVI scores include 
modeled results from just one of the 24 different 
models. When we excluded these modeled results 
from our CCVI scores, some of the CCVI scores 
moved up or down a score level (e.g., from highly 
vulnerable to moderately vulnerable or increase 
likely to presumed stable; see Table 4 and Fig. 4 
for an example). In other words, the results of the 
single model used in the CCVI have a large 
impact on the final CCVI score, which is 
worrisome since modeled results are extremely 


dependent on model algorithm and the choice of 
predictor variables. None of our species were 
scored as “increase likely” when we excluded the 
modeled results, despite that approximately half 
of the factors could be scored as decreasing 
vulnerability. Thus, the potential positive effects 
of climate change as measured by the CCVI were 
driven solely by increases in future range size 
predicted by SDMs. However, SDM results may 
be relatively reliable for predicting where suitable 
habitat will be in the future, compared to 
predicting declining suitability (Schwartz 2012). 
While recent advances in modeling techniques 
(i.e., consensus modeling, modeling clades instead 
of species, adding process models, factoring out 
spatial autocorrelation in occurrence data) may 
help create a model that better reflects the niches 
of plants, we expect that a quantitative modeling 
approach, in isolation, will never be sufficient to 
accurately predict the fate of rare species in the 
face of climate change. 

Topographic Complexity Analysis 

Our index of topographic complexity may help 
create a finer understanding of climate change 
vulnerability for our species. For example, 41 
species are ranked as highly vulnerable, but they 
have a range of topographic complexity scores. 
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Fig. 4. Maps of the predicted current (left column) and future (right column) ranges for two vulnerable species. 
For Piperia yadonii, the predicted range loss causes the overall CCVI score to increase, from highly vulnerable to 
extremely vulnerable. For Monolopia congdonii , the predicted range gain causes the CCVI score to decrease, from 
extremely vulnerable to moderately vulnerable. 


Thus, we can single out the highly vulnerable 
species that also have low topographic complex¬ 
ity scores as being especially vulnerable, especial¬ 
ly if they have low dispersal ability, with the 
mechanistic expectation that the local topogra¬ 
phy will not be sufficient to buffer a species from 


region-wide climate change by providing local 
refugia. An additional benefit is that elevation 
data is available at higher resolution than climate 
variables, allowing for finer suitability predic¬ 
tions. Furthermore, topographically homoge¬ 
neous places have been predicted to have faster 
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velocities of climate change, at least when 
compared at the worldwide scale (Loarie et al. 
2009). That said, the relationship between local 
landscape features and climate is complex and is 
just beginning to be described (Dobrowski 2011). 
Moreover, the interactions among topography, 
soils, soil water capacity, and microclimate on 
plant performance remains poorly described, 
despite the fact that our observations suggest 
plants are extremely sensitive to such interac¬ 
tions, at least in California and, more generally, 
in the mid-latitudes. 

CCVI for Rare Plants 

NatureServe’s CCVI is an excellent structure 
and transparent clearinghouse for information 
regarding climate change vulnerability. To our 
knowledge, it considers the most exhaustive list of 
extrinsic and intrinsic factors that may influence 
vulnerability, and also allows input of model- 
based results. Further, it is in use by many 
different groups, allowing for comparison of 
results. However, we’ve identified several prob¬ 
lems with the CCVI as it applies to plants. 

We were able to process only one species per 
eight-hour workday, a rate far too slow and 
expensive for most agencies to use for all the 
listed species in California. Our approach of 
subsetting a larger list based on rarity type had 
the advantage of possibly identifying particular 
combinations of range size, population size, and 
habitat specialism that cross-walk to climate 
change vulnerability, thus saving us the need to 
complete the CCVI for the remaining species. 
However, we found rarity type had no predictive 
power for the CCVI scores, and thus a detailed 
species-level analysis seems necessary to rank 
species with regards to climate change vulnera¬ 
bility. An alternative or complementary approach 
is to conduct other types of assessments, includ¬ 
ing vulnerability of landscape features (topogra¬ 
phy and connectivity) and habitats. These assess¬ 
ments can be completed relatively easily, and the 
results are perhaps more reliable, given that the 
connection of vulnerability scores to landscape 
features and habitats is less tenuous than the 
connection of vulnerability scores to species 
distributions and species ecologies. For example, 
most conservationists might agree that a well- 
connected landscape is less vulnerable to climate 
change than a fragmented landscape, but fewer 
might agree that a species with three pollinators is 
less vulnerable to climate change than a species 
with only one pollinator. 

Some attributes that are important to plant 
vulnerability are missing, including mating sys¬ 
tem (selfing vs. out-crossing) and pollinator 
specificity and efficiency. We recommend that 
different “flavors” of the CCVI be released in the 
future, at least one for animals and one for 


plants. Also, it is nearly impossible to complete 
the scoring for a given plant species, because 
information is simply lacking. When information 
is lacking, the guidelines often recommend 
scoring the species as neutral, while “unknown” 
is often the more appropriate score. Also, some 
of the scoring guidelines are too simplistic. For 
example, soil endemics are scored as more 
vulnerable to climate change than soil generalists, 
while this remains an under-addressed research 
question (Damschen et al. 2012). A related issue 
is that soil specificity should be assessed as a 
natural barrier; currently, it is only considered a 
life history trait. Finally, while the CCVI 
accounts for species interactions in a general 
sense, it does not explicitly take into account 
invasive species, which have major impacts on 
California plant diversity. Invasive species can 
become more virulent or less virulent depending 
on temperature and precipitation changes, and 
can greatly affect a species’ native habitat. For 
example, a recent study showed that climate 
shifts could increase the dominance of exotic 
species (Sandel and Dangremond 2012). 

Conclusions 

The information produced in our vulnerability 
assessments will be useful in identifying the most 
vulnerable rare plant species to climate change, 
which can then be carefully monitored. More¬ 
over, vulnerability assessments are an excellent 
way to identify knowledge gaps and to form new 
hypotheses about species distributions and cli¬ 
matic tolerances. Viewing multiple sources of 
information together, including the CCVI vul¬ 
nerability score, rarity ranking, topographic 
complexity, and a range of SDM results may 
give a broader picture of the overall vulnerability 
of a rare plant species to climate change. 
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